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ABSTRACT 


THE GALACTIC DISTRIBUTION OP CARBON MONOXIDE: 
AN OUT-OF-PLANE SURVEY 


Richard Stuart Cohen 


Galactic CO line emission at 115 GHz has been surveyed 
in the region 15°<S,<60® and -l°.5<b<l°.5 using the Columbia 
4-foot telescope (HPBW 8 arc minutes) . 

This survey, which comprises the first systematic out- 
of-plane data on CO emission from the Galaxy, confirms the 
finding of previous in-plane studies that CO is concentrated 
in a ring 6 kpc in radius. It provides the first determina- 
tion of the thickness of this molecular ring as a function 
of galactic radius and shows that CO is displaced from the 
conventional galactic plane. These results were arrived at 
by least-squares fitting the survey data to a circularly 
symmetric model of the Galaxy. The average half- thickness 
at half maximum of the molecular ring is 59 pc and- the aver- 
age displacement of the ring with respect to the b=0° plane 
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A comparison of the CO and HI distrihutions shows that 
there are marked differences in the distributions of these 
species r both radially and out of the plane. The CO distri- 
bution is similar to those of HII regions, supernova remnants, 
and other young galactic objects. It is argued that CO is 
probably the best available tracer of this component of the 
Galaxy. 

A detailed discussion of the antenna characteristics, in- 
cluding the radiation pattern and pointing characteristics 
is presented. 
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I . INTRODUCTION 

Since the discovery of widespread interstellar mole- 
cules in the late 1960's, it has become clear that molecu- 
lar clouds are an exceedingly important constitutent of the 
interstellar gas- They comprise a large fraction, probably 
over half, of the interstellar mass and are the raw material 
from which stars evolve. So far, most of the effort in mole- 
cular radio astronomy has been directed toward the study of 
well-known ob 3 ects such as HII regions, stars, and nearby 
dark clouds. Recently investigation has begun into the 
distribution of molecular clouds on a galactic scale. The 
initial steps in this investigation are naturally to deter- 
mine several of the overall characteristics of the molecular 
matter: its distribution both as a function of galactic ra- 
dius and of distance from the galactic plane; and the na- 
ture of its fine scale structure, for example its dumpi- 
ness and velocity dispersion. 

The best tool currently available for studying these 
problems is line emission by carbon monoxide. Although H 2 is 
by far the most abundant molecule, it has no dipole moment 
and cannot be seen at radio frequencies. It can be seen from 
satellites in ultraviolet absorption, but these observations 
are generally limited by interstellar extinction to within 
a distance of about 2 kpc from the earth. On the other hand. 
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emission from the first rotational transition of CO at 115 GHz 
has been found to be easily detectable in virtually all 
massive molecular ob;]ects. Although this line is usually 
found in clouds in which it is optically thick, these clouds 
seem to be distrubuted sparsely enough that there is no 
serious problem of obscuration of distant clouds by fore- 
ground matter. Thus, using CO as a probe, molecular matter 
can be detected throughout the Galaxy. 

This thesis reports the results of the first systematic 
out-of-plane survey of CO in the galactic disk. We will 
obtain the distribution of CO — and by implication the 
distribution of molecular matter generally — both as a 

i 

function of galactic radius and, in an important new result, 
as a function of distance from the galactic plane. 

Radial distribution of molecules 

Until now, the distribution of molecular matter had 
been systematically investigated only in the galactic 

plane. Two independent groups have made surveys of CO in 
the plane (Scoville and Solomon 1975; and Burton et al . 1975, 
Burton and Gordon 1976) and both have arrived at similar 
conclusions: excluding the galactic center region (R<600 pc), 

where intense emission is seen, the bulk of the CO emission 
is concentrated in a ring about 2 kpc wide and centered 
at about R = 5.5 kpc from the galactic center. This CO 
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distribution is similar to those of HII regions (Burton 
et al , , 1975), diffuse HII gas (Gordon and Cato, 1972), 
supernova remnants (Ilovaisky and Lequeux, 1972; 

Kodaira, 1974), and pulsars (Seiradakis, 1977), i.e. all 
the Population I objects whose galactic distributions are 
known with the exception of HI. The HI distribution is 
however very different; HI emission rises smoothly from 
near the galactic center to a peak at R = 14 kpc. 

There is also a striking correlation between CO 
emission and y-ray emission (Puget and Stecker, 1974; 

Paul ^ , 1976) . Because high energy y-rays (> 100 MeV) 

are produced by the interaction of cosmic rays with the 
nuclei of interstellar gas (Stecker, 1977) , the y-ray 
emissivity is expected to be proportional to the product 
of the interstellar mass density and the cosmic ray den- 
sity. Thus the correlation of y-ray emission with molec- 
ular emission, rather than with 21 cm emission, suggests 
that H 2 is the dominant component of the interstellar mass. 

It should be emphasized that before the present work 
the galactic CO measurements were made at b = 0 ° , while 
the data for the other tracers is integrated across the 
disk of the Galaxy. A proper comparison therefore 
requires extensive out-of -plane CO data. 
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Molecular Mass of the Galaxy 

To estimate the total molecular mass of the Galaxy from the 
in-plane CO surveys, it has been assumed that the average 
ratio of H 2 mass to CO luminosity for molecular clouds 
is the same throughout the Galaxy. In addition, as the 
CO surveys have been confined almost exclusively to the 
galactic plane, a thickness of 120 pc — estimated from 
a few out-of-plane measurements — was used for the entire 
galactic disk. 

Several methods have been used to determine the mass- 
to— luminosity ratio for molecular clouds. With a simple 
radiative transfer model for molecular clouds. Burton and 
Gordon (1976) used measurements of the intensities of 
the J = 1 -»■ 0 transitions of both ^^CO and 
to calculate the column density of CO. 

Then, taking the value of [C0]/[Il2] from local measurements 
to be constant over the entire Galaxy, they obtained the 
H 2 column density. Alternatively, the actual ratio of H 2 
column density to CO intensity in the direction of the 
galactic center can be measured (Stacker ^ 1975) , 

The total hydrogen column density — molecular and 
atomic — can be obtained from either infra-red extinc- 
tion or X-ray absorption data; the atomic component, as 
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measured at 21 cm, must then be subtracted to get the 
column density. Each of these methods gives 
about the same total galactic molecular mass: 

Q 

2x10 M_ . Unfortunately, they all involve 

so many uncertain assumptions that they may not give more 
than order-of-magnitude estimates. 

In principle, the mass can also be obtained directly 
from the data. The y-XELy emission is proportional 

to the cosmic ray density times the total gas density. 

Thus if a model cosmic ray distribution is assumed and 
if the distribution of atomic matter is known, then the 
molecular mass can be computed. Models with moderate 
(less than a factor of two) radial variation in the cosmic 
ray density give masses consistent with those derived 
from CO data (Stecker et , 1975) . However, the cosmic 
ray density is poorly known at large distances. If the 
hydrogen density could be more accurately determined from 
the radio data alone, the y-ray data could be used more 
effectively as a measure of the cosmic ray density in the 
interiors of molecular clouds. 

Summarizing, the best current value for the mass 

9 

of the Galaxy is 2 x 10 M^; however, this is really 

little more than an order-of-raagnitude estimate, T'^hen 
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detailed studies of distant molecular clouds have been 
made, it will be possible to better determine their 
molecular content and improve the accuracy of the mass 
estimate. 

Thickness of the Galactic Plane 

As we have seen, to determine the distribution it is 
essential to know the thickness of the molecular layer. 

In addition, the thickness is interesting in its own right 
as a probe of galactic structure. Investigations of 
the thickness of the stellar component of the Galaxy date 
from the beginnings of Galactic astronomy (see, e.g, Obrt, 1932). 
Two major results of these studies have been of great impor- 
tance. First, from measurements of the density and velocity 
dispersion of a single type of star as functions of height 
above the galactic plane, z, it has been possible to deter- 
mine the gravitational field and thus the total mass den- 
sity, both as functions of z (Oort, 1965; Woolley, 1965); 
and second, information on the scale heights of the various 
stellar types has been crucial in classifying the stars into 
populations (Blaauw, 1965) . 

In spite of the great success of these stellar pro- 
grams, they do suffer from a severe limitation: optical 

observations can generally be made only within about two 
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kiloparsecs from the sun. This restriction can be removed 
if radio data is used. We will briefly summarize the radio 
data available on the non-local thickness of the galactic 
disk. The radio observations can be divided into two 
broad categories: diffuse gas and discrete objects. The 

HII gas can be detected by hydrogen recombination line 
observations. (Actually there is some question as to 
whether these lines really come from diffuse gas or simply 
from weak HII regions.) Lockman (1976) and Gordon et al . 
(1972) have measured the half thickness at half maximum (z, ) 
of the HII gas at between 40 and 80 pc. However, these 
measurements are extremely time consuming and this result 
is based on a total of only 10 spectra. The distribution 
of HI gas, easily observed at 21 cm, has been thoroughly 
studied, but far from the sun it is very difficult to 
distinguish the diffuse HI gas from the HI clouds. Using 
a rather complex stochastic model of HI clouds distributed 
in a diffuse medium. Baker and Burton (1975) obtained 
Zj^ = 140 pc for the diffuse gas. In addition to the 
observational problems, there are dxfficult problems in the 
interpretation of the thickness data because the dynamics 
of the diffuse gas are strongly influenced by thermal, 
cosmic-ray and magnetic pressures. 


The effects of these pressures are reduced or eliminated 
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in the case of discrete objects. The discrete radio ob- 
jects that can be seen at large distances i-ncl-ude cold H-I 
clouds, HII regions, supernova remnants (SNR’s), pulsars, 
and molecular clouds. There is an enormous amount of HI 
data, but, as we have seen, to distinguish the clouds from 
the diffuse medium is extremely difficult. For R < 10 kpc, 
the Baker and Burton (1975) stochastic model gives an aver- 
age value of z, = 95 pc but no information on the variation 
of with R, 

Pulsars, SNR's and HII regions can be observed through- 
out the Galaxy, but the niaraber of objects bright enough to 
be seen up to 10 kpc away is less than 100 in each case. 

This is large enough to measure the average scale height 
of the Galaxy, but not adequate to measure the scale height 
as a function of galactic radius. For pulsars, z, = 330 pc 
(Seiradakis, 1977), much larger than for other population I 
objects. The large thickness is presumably the result of 
large velocities given to the pulsars at their formation. 

For SNR's (Ilovaisky and Lequeux, 1972) and HII regions 
(Burton et ad. , 1975) z-^ ^ 67 and 4 0 pc, respectively, consistent 
with other extreme population I objects. Ilovaisky and Leqeux 

did make a rough determination of the radial dependence of 

2 

the thickness of SNR' s and found that < z > ^ increases from 4 6 
pc for 4<R<8< kpc to 89 pc for 8<R<10 kpc. If the thickness of 
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either HII regions or SNR's is to be determined in more de- 
tail, many new objects will have to be detected. It is not 
likely that this will be possible in the near future. 

In contrast to the other radio objects we have considered, 
molecular clouds are extremely numerous. There are at least 
10,000 according to the cloud density estimate of Scovxlle 
and Solomon (1975) , and virtually all of these can be ob- 
served with existing telescopes. The molecular clouds are 
probably sufficiently heavy that, like stars, their motion 
is determined by gravity alone and not by the pressure of 
the interstellar medium (Mouschovias , 1975) . If the thick- 
ness and velocity dispersion of the molecular clouds can be 
deteimilned throughout the Galaxy, it may be possible to use 
them as a probe of the gravitational field on a galactic 
scale, just as stars have been used locally. 

Prior to the present work, the best published measurement 
of the thickness of the CO layer had been made by Burton and 
Gordon (1976) based on a single series of spectra taken at 
Si = 21° between b = -2° and +2°. By generating a number of 
synthetic velocity vs latitude contour maps of CO emission 
based on models with various thicknesses of the CO layer and 
then choosing the map that best matched the actual observa- 
tions, they estimated that the data were best fit with a 
Gaussian dispersion about the galactic plane of 50 pc. They 



10 


made no attempt to determine the variation of the thickness 
with ga-l-ae-feic rad-ius. 

There are several reasons for a more detailed 'deter- 
mination of the thickness of the molecular layer. First, 
the molecules are a major component of the interstellar me- 
dium, and a thorough investigation of their distribution is 
essential. Second, because of their large numbers and ease 
of detection, molecular clouds are the best tracer of the 
extreme population I component of the Galaxy. And third, 
because they move as pressure free objects, the molecular clouds 
may serve as a probe of the galactic gravitational field. 

Present Work 

To make an extensive survey of the molecular content 
of the Galaxy, the construction of a 4— foot radio telescope 
on the roof of the Pupin Physics Laboratory at Columbia Uni- 
versity in New York City was undertaken in 1972. Because of 
the short wavelength at which it is operated — 2.6 mm, the 
first rotational transrtxon of CO — the telescope, while 
quite small, has approximately the same resolution as the 
best that is available in any large scale 21 cm survey. The 
telescope was completed in December 1974, and a systematic 
survey of CO in the galactic disk was begun in November 1975 . 
The results of the first phase of that survey are the sub- 
ject of this dissertation. 
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II. APPARATUS 

The observations were made using the 4-foot telescope 
at Columbia University in New York City. The telescope 
system, diagrammed in Figure II-l, consisted of five main, 
elements: 1) the antenna, 2) the mount, 3) the front end 
which amplified the incoming microwave signal and converted 
it to a lower frequency, 4) the spectrometer or backend, 
and 5) the computer system which controlled the telescope 
and processed the data. Each of these elements will be 
described briefly below. More detailed descriptions can 
be found in the theses of those principally responsible 
for each element: details of the antenna characteristics 

and the mount alignment are in appendices to this thesis, 
the front end is described by H.I. Cong (1977), and the 
backend and computer systems are described by G. Chin (1977). 

The antenna, built by Philco Ford and located on the 
roof of the Pupin Physics Laboratories, is a Cassegrain 
reflector with an 8 arc min beamwidth (FWHM) at the CO 
frequency. It consists of a 4-foot (1.2 meter) main 
reflector with f/D = 0.375 and a 6 inch (15 cm) secondary, 
giving an effective system f/D = 2.8. The antenna was 
specified and was mechanically measured by the manufac- 
turer to have a surface accuracy of better than X/75 at 
the 2.6 mm CO wavelength, assuring that it would be 
diffraction limited. 
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Figure II-l, Telescope block diagram. 
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To verify that these specifications were actually met 
and that the antenna was accurately aligned after delivery, 
the antenna pattern was measured at the CO frequency 
of 115 GHz using a signal transmitted from a distance of 
420 meters. For comparison, a theoretical antenna pattern 
was calculated using scalar Kirchhoff diffraction theory. 
When the effects of the aperture blockage by the secondary 
mirror and its support structure were included, there 
was excellent agreement between the measured and calculated 
patterns (Figs. II-2 and II-3) . These results were 
adequate to demonstrate that the dish did in fact meet its 
manufacturing specifications - 

The effects of errors in feed horn positioning were 
also computed. It was found that the horn must be positioned 
to within 2 cm of the focus along the telescope axis and 
to within 11 cm transverse to the axis, and that the horn 
and telescope axes must be colimated to within 1°, These 
calculations were checked experimentally by deliberate 
maladjustment of the horn. In practice, the specifica- 
tions were easily met. A more detailed discussion of all 
the calculations and measurements is given in appendix B. 

The antenna was mounted on an altitude-azimuth drive 
built by the Physics department machine shop at Columbia 




Figure IX-2. Antenna pattern in horizontal plane with 
13 dS edge illumination. Solid lines are the calculated 
patterns; dotted lines are the measurements. 





Relative Power {db) 



Figure II-3, Antenna pattern in vertical plane with 
13 dB edge illumination. Solid lines are the calculated 
patterns; dotted lines are the measurements. 
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University- The position of the telescope was sensed by 
the computer system by means of 15 bit optical shaft 
encoders mounted on the telescope axes. Twenty times per 
second the computer read the encoder values and compared 
them with the desired values. If a correcetion was needed, 
the computer sent pulses to the stepping motors that con- 
trolled the telescope motion. 

Because the sun and moon are the only astronomical 
sources whose positrons can be readily measured at the 
CO frequency with a 4-foot telescope, the accuracy of the 
drive system could not be established over the entire sky 
by simple radio checks. Therefore, an optical telescope 
mounted on the radio dish and collimated with the radio 
axis of the telescope was used to check the pointing. 

The collimation was accomplished to within 0.25 arc minute 
by tracking the center of the sun with the radio axis and 
adjusting the optical telescope until its crosshairs were 
centered on the sun. Then, using the optical telescope, 
the pointing accuracy was measured against the positions 
of bright stars. 

In the initial star-pointing tests, corrections were 
included for atmospheric refraction, 'deviation of the 
azimuth axis from the true vertical, non-perpendicularity 
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o£ the altitude and aximuth axes, and non-perpendicularity 
o£ the optical and altitude axes. The tests established 
that the two non-perpendicularity corrections were small 
enough they they could be dropped with no effect on the 
pointing accuracy. The magnitudes of the other corrections 
were determined from preliminary tests, and then a final 
check was made. 

Figures II— 4 and II-5 are plots of the residual 
errors for the 31 stars included in the final test. As 
these graphs show, it was possible to point the telescope 
to within 2 arc minutes of any desired position. The 
RMS error was 1.1 arc minutes. We estimate that the error 
was primarily due to three effects of about equal mag- 
nitudes: 1) quantization errors in the encoders and the 
stepping motor drive, 2) xrregularities in the encoders, 
and 3) lack of rigidity in the mount. While the 2 arc 
minute (k beamwith) accuracy is adequate for present 
purposes, significant improvements are planned for the 
near future. 

The pointing accuracy was confirmed by direct radio checks 
against the position of the sun that were made several 
times a week during the course of the observations. The 
accuracy of the telescope was further corroborated when 
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Figure II-5. Residual pointing errors; azimuth ajcis. 
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many sources mapped with it were later mapped with higher 
resolution telescopes. No discrepancies in position were 
ever found. Additional details of the telescope alignment 
procedure and the pointing tests are given in appendix A. 

The front end was a double sideband superheterodyne 
receiver constructed by H.I. Cong. The incoming radiation 
was collected by a feed horn, and using a resonant ring 
injection cavity and a Schottky barrier mixer built by 
A.R. Kerr, was mixed with the signal from a phase locked 
klystron local oscillator to produce a 1.4 GHz intermediate 
frequency signal. This signal was amplified by a Micro- 
mega room temperature parametric amplifier with a 50 K 
noise temperature, further amplified by a low noise tran- 
sistor amplifier, and then down converted to the second 
intermediate frequency of 150 MHz. (Por details of the 
front end see Cong, 1977.) The complete system had a 
measured single sideband noise temperature of 1400 K 
referenced to the feed horn. The 150 MHz signal was 
amplified and sent to the spectrometer located in the con- 
trol room on the floor below. 

The spectrometer contained 40 three stage Butterworth 
filters, each with 1 MHz resolution full width at the 
half power points (corresponding to a Doppler velocity 
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resolution of 2.6 km sec”^ at the CO frequency), giving 

_1 

a total bandv?idth of 4 0 MHz (104 km sec ) . The outputs 
of the filters were detected and sent to analog integrators. 
Every 50 msec the outputs of the integrators were digitized 
by a 12 bit analog- to-digital converter and sent to the 
computer for further processing. 

The computer system had two main functions: 1) to 
control the position of the antenna, and 2) to collect, 
process and store the data from the spectrometer. It 
consisted of a Data General Nova 1220 computer, a fast 
fixed head disk, 2 IBM compatible tape recorders, a paper 
tape reader and punch, a Tektronix graphic display termi- 
nal and hard copy unit, and a clock unit that provided* 
time-of-day and synchronization signals. The programming 
was done by G. Chin and G. Tomasevich with assistance from 
the author. Details of the programs and hardware are 
given in Chin (1977) . 
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III. OBSERVATIONS 

Observations of the Galaxy were made every 2°. 5 in 
galactic longitude from 15° to 45° (except 22°. 5) and every 
5° from 45° to 60°. At each longitude, observations were 
spaced 0°.25 (about two beamwidths) in latitude. In most 
cases the survey was extended above and below the plane 
until two successive positions showed no signs of CO 
emission. 

The previous in-plane surveys have found almost no 
CO emission at negative radial velocities in the longitude 
range of this survey. Because our spectrometer was not 
wide enough to make observations of both positive and neg- 
ative velocities simultaneously, it was therefore decided 
that, for now, the previous null results would be accepted 
at negative velocities, and observations were thus confined 
to positive velocities. For I ^ 40° neither HI nor CO 
emission has been observed at velocities with respect to 
the local standard of rest (v ) greater than 100 km sec 

IjDXv 

At these longitudes it was therefore possible to include in 
the 104 km sec ^ bandpass of the spectrometer the entire 
range of velocities at which emission was expected. Inside 
.£_=:__40- --the -receiver was set to i-ncl-ude the highest velocity 
at which significant emission was found by Burton et al. (1975) . 
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At these latitudes the low velocity (generally <10 km sec ) 

emission, presumably originating in local CO, lay outside 
the range of the spectrometer, and was consequently not 
observed. After the completion of this survey, Gordon and 
Burton (1976) published a high sensitivity in-plane survey 
that showed some emission at velocities higher than those 
covered in the present survey at Z = 15°, 17°. 5, and 20°. 

The omission of this high velocity emission will only affect 
the conclusions concerning the inner 4 kpc of the Galaxy. 

All the spectra were obtained by position switching, 
i.e., by observing during alternate 30 sec periods at the 
source position and at a comparison position. The compar- 
ison positions ("offs”) were checked with considerable care 
to assure that they were free of any CO emission stronger 
than 0.3 K (less than half the 3o noise level of the survey 
itself) . This was done initially by the frequency switching 
technique in which a comparison spectrum is generated by 
shifting the frequency of the receiver by half the bandwidth 
of the spectrometer instead of by moving the telescope to 
an emission free position. While this scheme usually produces 
baselines that are inferior to those obtained by position 
switching, it is ideal for finding emission free regions 
because it does not involve a second, possibly "contaminated", 
off position. The comparison positions were further examined 
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by position switching against each other. The final com- 
parison positions, listed in Table III-l, were typically 
5“ above and below the galactic plane and 5° apart in longi- 
tude. 


The receiver was calibrated against a room temperature 
blackbody by a chopper wheel technique similar to that 
described by Davis and Vanden Bout (1973) . Their method, 
which corrects for the beam efficiency and for the atmos- 
pheric attenuation in a single layer atmosphere, was refined 
by using a two layer model; an upper layer of ©2 and a lower 
layer of H 2 O. The optical depth of the O 2 was calculated 
from a standard atmosphere to be 0.21 at the zenith, while 
the optical depth of water was derived from antenna tippings 
made at least once every six hours. (See Chin, 1977, for 
details of the calibration techinque.) Observations were 
generally made on cold, dry winter days when the optical 
depth of water was between 0.15 and 0.25 at the zenith. 

For a given source this calibration procedure gave 
corrected antenna temperatures, T*, that were independent 
of elevation down to the lowest angles above the horizon 
at which observations were practical (about 20*'). Figure III-l 
sho ws a series of measuremL^nts , made at different altitudes, 
of the calibrated intensity from the position of the most 
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TABLE III-l 


Emissionr-free comparison posxtions 


1 r) 

15 

15 

20 

25 

30 

35 

35 

40 

40 

43 

45 

45 

50 

50 

55 

55 

60 

60 


b 

-5 
+5 
+5 
-5 
.-5 
-5 
+5 
-5 
+5 
-2 
+ 2 
+3 
-3 
+2 
-5 
+ 5 
-5 
+5 
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intense emission in Orion A, integrated over the line and 
normalized to the value measured at transit. It can be seen 
that over the altitudes at which observations are normally 
made, the calibration procedure compensates extremely well for 
atmospheric attenuation. 

The absolute calibration was checked .on several astro- 
nomical sources against data taken with the NRAO 36-foot 
telescope at Kitt Peak. The Kitt Peak results were spatially 
smoothed to the resolution of the 4-foot telescope and the 
intensities, integrated over linewidth, were compared. Agree- 
ment was found to be within 20% - as good as can generally 
be obtained by current techniques. 

Because the atmosphere contributes different amounts of 
noise at the signal and comparison positions, the spectra were 
offset from zero by up to 30° K. In most cases the offset was 
reduced to less than 10° K by averaging two scans whose com- 
parison posxtions were displaced roughly equal amounts in alti- 
tude on either side of the sxgnal position. In all cases it 
was necessary to subtract from the spectrum a straxght line 
that was least- squares fit to parts of the spectrum free of 
emission (usually the ends) . Tests on sky regions free of 
any CO emission demonstrated that this method removed all 
baseline effects to well below the noise level of this survey. 
(See, for example, figure nx-3^ 2- - 40°, b = 3°. 25). 
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Figure II 1-2, a summary of the observational data, shows 
the corrected CO antenna temperature., T*,. a.t various gala.^ctic 
longitudes as a function of galactic latitude and radial 
velocity in the local standard of rest. This figure repre- 
sents a total of 179 spectra. Each spectrum is the result of 
between 20 and 40 minutes of integration time. The integra- 
tion time was adjusted, depending on atmospheric conditions, 
to give a 3a noise level of about 0.7 K. A sample of 
individual spectra are shown in Figure III-3. 

Because we have not mapped closely in latitude, it 
has not generally been possible to demonstrate clear rela- 
tionships between the CO clouds shown in Figure III- 2 
and known astronomical objects. However, we have made 
preliminary identifications of a few of the more prominent 
features, indicated by letters dn the contour maps and 
described below. 

& 15° ; Feature A is M17 . This HII region has a large 

associated molecular cloud that has been thoroughly mapped 
in CO by Lada (1976) . 

I = 17° , 5 : Feature B is probably related to Ml6, a 

large HII region at G17.0 +0.8 . Georgelin and Georgelin 
(1976) give the H109a velocity as 24.5 km sec ^ and the 
Ha (optical) v elocity as 28.4 km sec 
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Figure Distribution of CO emission in the galactic 

disk. The contour interval is 1 K antenna terrperature (cor- 
rected for beam efficiency and atmosphere) . The area of each 
plot shows the region actually surveyed in both velocity and 
latitude. The angular resolution indicated is that of the survey 
(2 beamwidths) . Two additional longitudes were surveyed but 
showed no enission at the first contour level; S, = 55° fran 
b = -0°.75 to b = +0°.75 and v^„^ = 0 to 104 km sec ; and 

J_oK 

£ = 60° fron b = -0°.5 to +0°.5 and v^„„ = 0 to 104 km sec 

IjoK 

The letters indicate features that are described in the text. 





CO 
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Jl= 20° : Feature C is the HII region Gl 9. 7-0.2 which has 

an Hl09a velocity of 43.4 km sec~^ (Georgelin and Georgelin, 
1976) . 

i = 25° ; I^eifenstein et. al . (1970) note that this is 

a complex region. They list two Hl09a sources near 
feature D, namely G24. 6-0.2 at 109 km sec and 
G24, 8+0.1 at 114 km sec They give a velocity of 

60 km sec“^ for W42 at G25.4-0.2; this source is probably 
related to feature E. Clark and Caswell (1976) show a 
SNR at G24.7-0.6 with a diameter of 15 arc minutes, large 
enough to suggest that it may be related to an object on 
this map. 

A = 27°. 5 : Georgelin and Georgelin (1976) give an Hl09a 

velocity of 98 km sec ^ for the HII region at G27.3-0.2, 
probably feature F. 

I - 30° : Feature G is related to W43, a complex region 

containing many sources. The CO is most likely associated 
with the HII region at G29. 9-0.0 at an H109a velocity of 
96 km sec ^ (Reifenstein et. al. , 1970) . There is also a 
SNR at G29. 7-0.2 (Clark and Caswell, 1976). 

^ = 35° : There are numerous sources in this region. 
Feature H is a local dark cloud that has been mapped in 
H^CO and OH by Myers (1975) . Feature J is the complex 
W48 region which is related to a large HII region. 

Feature K is near the SNR W44 located at G34.6-0.5. 



^ = 40*^ ; Feature L is near W45 at G4 0.5+2. 5 and is 


-probabl-y- the +I-FI region -S 7-6 {-Sharp-iess, 1-959-)-. S-7-5 i-s - 

at G40.2+1.5, and it is possible that a related molecular 
cloud appears on this map. 
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IV. DISTRIBUTION OF CO EMISSION 

Before proceeding with a close quantitative analysis, 
we will attempt to determine the major features of the 
galactic distribution of CO directly from a careful exam- 
ination of the contour maps (Fig, III-2 ), First, we 
see that there is very little emission for S, > 45® . 

Since observations in these directions only sample points 
more than 7 kpc from the galactic center, we can conclude 
that for R > 7 kpc the amount of CO is very small. 

Next, in the inner part of the Galaxy we find a 
discrepancy between the maximum velocities of the 
HI and CO emissions. For 15° < I < 25° the maximun 
HI velocity increases with decreasing longitude, whereas 
the maximum velocity at which significant CO emission is ob- 
served decreases; the high velocity CO is extremely weak or 
missing. Because of the differential galactic rotation, the 
highest velocity emission at a given longitude originates 
closest to the galactic center. Thus the missing emission 
for < 25° implies a lack of CO inside of the £= 25 ° tangent 
point, i.e. for R < 4 kpc. So, we can see that the CO 
emission is concentrated in a molecular ring between about 
4 and 7 kpc from the galactic center. 
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To estimate the out-of-plane dxstribution of 
CO emissxon, all the spectra at the same latitude were 
integrated over frequency and averaged together. The 
result (see Figure ly^l ) shows that the CO plane is 
1°.2 thick (FWHM) and is displaced 0°.3 towards negative 
latitude. Taking R = 6 kpc as the average radius of the 
molecular ring and ignoring, for now, radiation from the 
far side of the Galaxy, we find that the distance from the 
sun to the ring is between 4 kpc (toward the galactic 
center) and 8 kpc (toward the tangent point) . Using 6 kpc 
as typical of the distance to the ring, the angular thick- 
ness of the molecular disk corresponds to about 130 pc (FWHM) 
and the displacement to 30 pc. 

In order to analyze the survey data more quantitatively 
we will assme that although local, observations indicate 
that CO is generally found xn clouds in which it is optic- 
ally thick, the Galaxy is optically thin in the sense that 
the clouds have a sufficiently sparse distribution that 
radiation from distant clouds is not obscured. While there 
are usually several clouds along any given line of sight 
passing through the molecular ring, differential galactic ‘ 
rotation will generally Doppler shift the radiation from 
distant clouds to frequencies at which more nearby matter 
is completely transparent. The best evidence' that the 
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molecular clouds are in fact sparse enough not to obscure 
each other i_s the appearance of the high resolution spectra 
taken with the NEAO 36-foot telescope as a series of well 
defined features separated by velocities with no emission 
(Burton et ad. , 1975) . If there were many clouds along the 
line of sight with similar Doppler velocities, the emission 
would not be likely to fall to zero between clouds. 

To convert the Doppler velocities of the CO clouds into 

galactic positions, a simplified model of galactic rotation 

has been used. Specifically, it is assumed that the large 

scale structure of the Galaxy is azimuthally symmetric 

and that the only motion is circular rotation. Although the 

azimuthal symetry of the model will supress evidence for 

spiral arms, density waves, etc., it will have only a 

small effect on the radial distribution derived here. The 

fluctuations in the observed rotation curve, sometimes 

attributed to non-circular "streaming" motions, are generally 

—1 

less than 10 km sec . If the simplified circular rota- 
tion model is used, these non-circular motions can produce 
as much as a 2.5 kpc error along the line of sight; how- 
ever, the corresponding error in galactocentric distance 
is never more than 1.3 kpc over the region of sky covered 
in the pre s ent survey. For 4<R<8 kpc, the region where 
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almost all the CO emission is concentrated, the error is 
less than 0.7 kpc. 

Because of its analytic conveni'ence , the rotation curve 
from Burton (1971) has been adopted. This curve, derived' 
from 21-cm data, gives the velocity as a function 
of galactic radius to be 

V(R) = 250.0 +4.05 (R^ - R) - 1.62 (R^ - R) ^ (IV-1) 

—X 

where V(R) is in km sec , R is the galactic radius in kpc, 

and R =10 kpc is the distance of the sun from the galactic 
o 

center. Por 3<R<10, this curve differs from the Schmidt 

—1 

curve by less than 3 km sec and from the curve derived 
by Gordon and Burton (1976) from CO data by less than 
2 km sec These discrepancies are less than the obser- 

vational uncertainties in the rotation curve due to small 
scale fluctuations and the criterion used to determine the 
maximimi Doppler velocity at a given longitude. 

In the region interior to the solar circle covered by this 
survey/ it is impossible, even with simplifying assumptions 
about galactic rotation, to assign a unique position in the 
Galaxy to each observed feature; along a given line of sight 
the two points with the same distance from the galactic 
center have the same velocity component towards the sun. 

The galactocentric distance of a cloud is thus uniquely 
determined, but the distance from the sun is ambiguous. 
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Since it is impossible to establish the locations o£ 
individual clouds, we have instead fit the data to simple 
stochastic models. The CO is assumed to be contained in 
many small, randomly distributed clouds. We take the 
average number of clouds in the volume element dV, rad- 
iating in the J = l->0 line of CO with an intensity between 
B and B+dB in the frequency interval between v and v+dv , 
measured in the sun’s local standard of rest, to be given 
by a distribution function 

g(R,z,B,v) dV dB dv (IV-2) 

where R is galactic radius and z is height above the gal- 
actic plane. If the clouds are all of the same radius, a, 
and if g does not vary significantly over a distance 
comparable to a, then the average number of clouds inter- 
sected by a given line of sight between a distance s and 
s+ds from the sun will be given by 

f (R,z,B,v) ds dB dv = TT g(R,z,B,v) dV dB dv (IV*-3) 

where R and z are now functions of s and the galactic 
coordinates Z and b. 

If we ignore the distant wings of the line and assume 
that the CQ clouds ar^_QpJ:i_c.al.lyL_ thick over the entire CO 
line, then we will see only the nearest cloud along the 
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line of sight emitting at a given frequency. Therefore, 
the expectation value of the observed intensity will be 

M 00 

I()l,b,v) = / Q Iprobability there is no cloud nearer 

than s to the sun] x {/^ [probability 
of cloud of brightness B at a 
distance s] B dB } ds 
= [ exp(-T^(s)) ] X e^(R,z) ds (IV-4) 

where t^(s) = f{R,z,B,v) dB ds is the average num- 

ber of clouds out to a distance s, and £^(R,z) = 

B f(R,z,B,v) dB . Equation IV-4 is the same as the 
solution of the equation of transfer in local thermody- 
namic equilibrium with an effective optical depth and 
an effective volume emissivity . 

If, as we have assumed, there is a low probability 
for more than one cloud to be along a given line of sight 
at a given frequency, i.e. t << i, then we can ignore the 
exponential term in Equation IV-4 and treat the Galaxy 
as effectively optically thin. Then we have 

I(il,b,v) = /“ e^(R,z) ds . (IV-5) 

Finally, we define the effective integrated volume emis- 
sivity as 
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e(R,z) = / e^(R,z) dv , (IV-6) 

where the integral extends over the GO line. It is this 
last function that we attempt to determine. 

Because it is usual in radio astronomy to express 
intensity in temperature units as T* = ( ) I , 

we will use compatible units for the effective emissivity 
and define 

e* = (c^ / 2kv® ) e (*K Mhz kpc . (IV-7) 

The effective emissivity will always be expressed in these 
units in this thesis and we will generally drop the 
asterisk. 


For our first stochastic model of the Galaxy, the 
effective integrated volume emissivity was assumed to be 
a gaussian in z, specifically. 


e (z ,R) 


Sq exp 


- (z-Zq) 2 (In 2) 



(IV- 8) 


where z is the distance above the galactic plane, R is the 
distance from the galactic center, and the displacement z^, 
the half-thickness at half power z, , and the central 
emissivity e^, are all functions of R that characterize 
the CO distribution. 
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To obtain the most likely values of the parameters 
Zj^(R), and z^(R) the Galaxy was divided rnto concentric 
rings of width 0.5 kpc. Initial estimates of the three 
parameters were chosen for each ring. Then, using the model 
emissivity function (Eq. IV-8 ) and the assumption of optic- 
al thinness, model spectra were computed for each latitude 
and longitude that was actually observed. The model spectra 
were compared with the observed spectra and the parameters 
were varied in each ring to obtain the best least-squares 
match of the model and observed spectra. Details of the 
fitting procedure and the computation of model spectra are 
given in Appendices C and D respectively. In rings with 
little or no CO there was insufficient data to fit all the 
parameters accurately. In that case the CO emission was 
assumed to be symmetric about the galactic plane, i.e., 
was constrained to zero. The resulting parameters as 
functions of R are shown in Figure IV-2 (left) , 

Proper convergence of the fit was assured in several 
ways. First, the fitting, v;hich was done by an iterative 
method, was continued until the change in each of the par- 
ameters from one iteration to the next was less than 2% of 
the calculated error in the parameter. Second, the entire 
fitting was repeated several times using different initial 
estimates for the parameters, including = 0.2 and 1; 
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R(kpc) 

Figure IV-2. Results of fitting data to a gaussian 
in z. Left; Displacement, half-thickness, and cen- 
tral emissivity of the CO disk as functions of galac- 
tocentric distance. Right: Half-thickness and 
central emissrvity with the displacement constrained 


to zero 
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z, = 10 pc, 50 pc, 200 pc, and 1000 pc; and .z =0 pc and 
^ o 

1000 pc. In every case in which the fit converged, it 
converged to the same value. 

Finally, as an overall check on the algorithm, the 
fitting was tested against a simulated version of the 
Galaxy, made up of randomly distributed clouds. The 
probability density for an individual cloud to fall at the 
^galactocentric cylindrical coordinates (R,0,z) was given 
by 

1 / (2 ir ) exp( / 20 ^^^ ) 

X exp ( - (z-*z ) ^ / 2a ^ ) x U(0) , (IV-9) 

o z 

where U(0) is a normalized uniform density between 0® 
and 360°, = 6 kpc, = 40 pc, and a^= 60 pc (corres- 

ponding to z, = 85 pc) . The clouds were all taken to 
have the same temperature, size, and linewidth. The size 
and linewidth were taken to be sufficiently small compared 
to the resolution of the telescope that the clouds can be 
treated as essentially pointlike with delta function emis- 
sion lines. 

From the simulated galaxy, spectra were constructed 
corresponding to what would be observed with the velocity 
and angular resloutions of the 4-foot telescope. These 
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spectr.a., which were saved by the computer in exactly the 
same format as the actual spectra, were then input to the 
least-squares fitting program as if they were real spectra. 
Figure IV-3 shows the results of the fitting for a sim- 
ulated galaxy. In this case there v/ere typically 15 clouds 
per resolution element in the portions of the spectra 
that correspond to a galactic position near the center of 
the simulated molecular ring. (Models with a much smaller 
number of clouds are considered in Appendix C.) It can 
be seen that the fitued values match the input distri- 
bution quite well and that the calculated error bars are 
a good indication of the actual errors . 

The actual data were also fit to a second, more restric- 
ted, model in which was constrained to zero everywhere. 
For R < 5.5 kpc the resulting (R) and (R) , shown in 
Figure IV- 2 (right) , do not differ significantly from 
the results of the previous fit. However, at greater 
radii, where the galactic plane is depressed according to 
Figure IV-2 (left) , there is a slight increase in z. . 

This is expected since the CO is typically further from 

z = 0 than from the true z . 

o 
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R(kpc) 


Figure IV-3* Results of least-squares fitting to sim- 

ula'ted galaxy. The dashed lines { ) indicate the 

underlying distribution used in simulating the galaxy. 

The dash-dot lines show the original parameter 

estimates used in the fitting. 
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Several checks were made to verify that the offset in 

shown in Figure IV~2 was not the result of a single large 

object far from the plane. The fitting was repeated twice, 

in one case ignoring the data for A = 35®, the latitude 

* 

with the largest concentration of emission below the plane, 
and in the other case the data for A < 25®, where the 
depression in the plane is most apparent to the eye in 

the contour plots (Fig. III-2) . In neither case did the 

offsets (see Figure IV-4) differ significantly from those 
obtained when all the data were included. 


Finally, the model emissivity function was changed to 
an exponential: 


e(R, 2 ) = Eq exp 


[z- 2 q| (In 2) 
* - 

ZJ3 


(IV-10) 


where e^, z^, and are again functions of R. The best 
values of the parameters were determined as before with 
the results as shown in Figure lv-5. 

If we compare these results with those from the gaussian 
model, we find that there is very little change in the shape 
of the curves. — We..-will now make a more -quantati-ve- com- 
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R(kpc) 

Figure IV-5. Results of fitting the data to an 
exponential in z. Left: Displacement, half-thickness, 
and central emissivity of CO disk as functions of 
galactocentric distance. Right: Half -thickness and 
central emissivity with the displacement constrained 


to zero. 
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parison of the models with three parameters (Figs. IV-2 
left, and IV-5,left) in the region 4.5<R<7.5, where all 
the parameters can be determined in both the gauss ian 
and exponential models. There is no significant change in 
the offset between the models: for the gaussian model, the 
average offset is <z^ = -24.3 pc, while for the exponen- 
tial model we get <z > — -25.5 pc. There is however a 

significant change in the Zj^ curve. Computing the average 
value of the ratio of for the gaussian model to that for 
the exponential model, we obtain < ^x, / 2 , > =1.45. 

Similarly, for the central emissivity we get 

< e /a > = 0,84. If we define the surface brightness 
O/ g o , e - 

of the Galaxy as 


e(RrZ)dz , (IV-11) 

then for the gaussian model 

= ( 7 r/ln 2)^ z, e , (IV-12) 

g =50 

and for the exponential model 

<y = (2/ln 2) z, e . (IV-13) 

e ^ o 


The value of 


< a / a > = 
e g 


/ 0^ averaged over radius is 
0.88 . Thus the change in z, between the 


two models is just about consistent with identical surface 
brightness . 
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The error bare in Frgures IV-2, -3, -4, and -5 

indicate the range over which the parameters may be varied 
without significantly affecting the quality of the fit as 
measured by , the mean square deviation of the model 
spectra from the observed spectra. {See Appendix C for 
details.) If the model matched the data very well, and 
were thus small, then a small change in a parameter would 
result in a significant fractional change in and con- 
sequently the error bars would be small. Thus rhe error 
bars measure, in effect, the inability of the model spectra 
to exactly reproduce the observations. Deviations of the 
observed spectra from the model spectra result primarily 
from random fluctuations in the distribution of molecular 
clouds. In principle, noise and inadequate baseline removal 
could also contribute to the deviations, but actually, they 
are insignificant for R > 3.5 kpc. 

Summarizing this section, by fitting the survey data 
with simple models of the distribution of CO emission, we 
have determined the thickness, the offset from z=0,- and 
the central eraissivity of the CO disk of the Galaxy. The 
results are insensitive to the detailed form of the 'model 
emissivity function. They are also not strongly influenced 
by the exclusion of a small part of the data.- Lastly, we 
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note that the fitting has confirmed the qualxtative esti- 
mates made at the beginning of this section concerning 
the existence of a molecular ring between about 4 and 7 kpc 
from the galactic center, and the thickness and offset of 
the CO layer. 

In view of the similarity of the results obtained by 
using either an exponential or a gaussian function for the 
model z-distribution of CO emission, our analysis alone 
presents no compelling reason to adopt one form over the 
other. However, there are arguments favoring the gaussian 
distribution: first, it is predicted from simple hydro- 
dynamic theory (Oort, 1965) and second, it seems to be 
a good representation of the distribution of HI clouds 
(Schmidt, 1957) . We have chosen therefore to use the 
results of the gaussian model (Figure IV- 2 , left) in 
the succeeding discussion. 
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V. DISCUSSION 


The major results of the analysis of Chapter IV are 
the thickness of the galactic CO disk, its offset from 
the conventional galactic plane, and its central emissivity 
— all as functions of galactic radius. In this chapter 
we will discuss these results and, in particular, their 
relationship to previous observations of the Population I 
component of the Galaxy. 

Central emissivity 

The central emissivity, (R) , (see Fig. V-1 which 
repeats the results of Fig. IV-2 ) is strongly peaked in 
a molecular ring lying between 5 and 7 kpc from the galac- 
tic center. Although a similar result had been obtained 
by the in-plane surveys, our conclusion is considerably 
stronger in that it is clear that we have properly located 
the position of the CO emission with respect to b = 0° . 

When our survey is compared in detail with the in-plane 
surveys, we find that if the emissivities measured by 
Gordon and Burton (1976) and Scoville and Solomon (1975) 
are normalized to our results in the molecular ring, then 
the Eq(R) curves agree reasonably well for R > 4 kpc (see 
Pig. V-1 ) . However, inside this radius, where the 
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R(kpc) 


Figure V-l, Emissivity of CO vs galactic radius. 
Filled circles; present work. Open triangles: 
Gordon and Burton (1976) sealed down by a factor 
of 4,3 to agree with the present work in the mol- 
ecular ring (see discussion in text) , 
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emission is weak, there is considerable disagreement: the 

Gordon and Burton -results are about -double those -of - - - 
Scoville and Solomon, while our survey shows no emission. 
Although Scoville and Solomon used the Schmidt rotation 
curve; Gordon and Burton used a curve derived from CO data; 
and we have used the curve given in Eq. IV-1, we have 
shown in Chapter IV that these differences in the choice 
of rotation curve are unlikely to cause large errors in 
the derived radial distributions. As was noted in 
Chapter III, Gordon and Burton found several clouds at 
5 , ^ 20® at velocities higher than those covered in our 
survey. It is the emission from these clouds that appears 
to be the source of the discrepancy. 

While it is, of course, desirable to have as wide a 
velocity coverage as possible and ‘the Gordon and Burton 
survey is thus somewhat more complete than ours for R < 4 kpc , 
one can not conclude that the Gordon and Burton results for e 
are necessarily a more accurate representation of the 
actual CO distribution. If, as was assumed by all three 
surveys (Gordon and Burton, Scoville and Solomon, and this 
thesis) , the galactic rotation were purely circular, and 
the Galaxy were circularly symmetric, then similar results 
-would— be -obtained us-ing -data— from- any ve-lociti-es— and—l-on- 
gitudes covering a given range of galactocentric distances. 
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The fact that the surveys obtained different results by 
observing at velocities and longitudes that were different 
but were predicted by the rotation curve to correspond 
to the same galactic radius suggests that either the rota- 
tion curve does not actually predict the correct radius 
or that the Galaxy is not circularly syimnetric . In fact^ 
we know from 21 cm data (Burton, 1974) that for R < 4 kpc 
there is a serious breakdown in the circular rotation 
hypothesis and probably also in circular syiranetry- We 
therefore conclude that the simplified models used by all 
observers for analyzing CO emission data are simply inade- 
quate for R < 4 kpc and that all the the results in that 
region must be considered highly speculative. 

For R > 4 kpc, the situation is quite different. There, 
as was demonstrated in Chapter IV, the elimination of the 
data from a wide range of longitudes has only a nominal 
influence on the derived distribution. In addition, 
all the CO surveys agree reasonably well in that region in 
spite of their different samplings in longitude. The 
21 cm data also indicate that our model is reasonable for 
R > 4 kpc. As was discussed in Chapter IV, the small 
scale fluctuations in the observed 21 cm rotation curve of 
about 10 km sec , presumably caused by non-circular 
motions, are not important for our purposes. The well- 
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known deviation between the northern and southern hemi- 
sphere rotation curves,, also caused by non-circular 
motions, is again about 10 km sec ^ and, for present pur- 
poses, indicates only a minor deviation from circular 
rotation . 

Because very different beamwidths were employed, we 
can not expect the spectra from this survey, done with a 
4-foot telescope, to agree in' detail with the spectra from 
the previous in-plane surveys, all of which were done 
with the NRAO 36-foot telescope; a 4-foot spectrum should 
be approximately the average of all the 36— foot spectra 
taken within the 4-foot’s beam. If the 4-foot's resolution 
is small enough that at a given velocity there is no sys- 
tematic large-scale variation in the CO distribution over 
its beam, then a 36-foot spectrum taken at the center of 
the 4-foot's beam will be typical of all the spectra that 
cover the beam, and the antenna temperatures recorded at 
a given velocity and position on the 4- foot and 36-foot 

f 

telescopes should agree on the average. For b = 0° , 
antenna temperatures from the tv70 telescopes can in fact 
be compared since the 4-foot beam corresponds to only about 
50 pc even at the most distant region surveyed (about 

2 0 kpc away) — consi der ably less than the 120 pc (FWHM) 

width of the CO layer. When integrated over all velocities 
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to eliminate the effects of different velocity resolutions/ 
we can expect that on the average the 4— foot and 36— foot 
spectra will agree. 

When this comparison is actually made, no two surveys 
agree; our results fall roughly a factor of 1.4 below those 
of Scoville and Solomon (1975 ) , a factor of 2 below those 
of Burton et. al. (1975) , and a factor of 3 below those of 
Gordon and Burton (1976) . Since the absolute calibrations 
of the 4-foot and NRAO telescopes agree reasonably well 
(see Chapter II) , and all the surveys done at NRAO used 
the same calibration procedure, this discrepancy does 
not appear to be an effect of calibration. 

The most likely explanation is the existence of an 
apparent weak substratum of CO emission covering much of 
the positive velocity region permitted by the rotation 
curve, and contributing most of the integrated intensity 
in the Gordon and Burton high sensitivity survey, but 
only partially detected in the other surveys. This emission 
has an antenna temperature of 1.2 K and, for S, < 35°, 
covers about 80 to 100 km sec ^ of each spectrum, con- 
tributing about 100°K,kra sec ^ to the integrated antenna 
temperature — about 60% of the average total integrated 
antenna temperature. Subtracting this substratum from 
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Gordon aiid Burton's results reconciles them with ours to 
within about 20% . 

The reason that the weak emission was not fully detected 

in all the surveys is probably related to their differing 

velocity coverages: the present survey covers 104 km sec ^ 

in each spectrum, Scoville and Solomon covered 166 km sec 

and Burton et. al. and Gordon and Burton covered 333 km 
-1 

sec . In both the present survey and in Scoville and 
Solomon, because the weak emission covered such a large 
fraction of the spectrum and had an intensity roughly -at 
the detection limit, it was presumably indistinguishable 
from the instrumental baseline and was subtracted out 
in the data analysis. ' Although we have not detected this 
weak emission, and the present survey is hence biased 
toward the more intense molecular sources, the close 
agreement in the radial distribution of all the surveys 
strongly implies that the clouds we are seeing do provide 
a reasonably faithful measure of the overall CO distribution. 

Thickness of the CO layer 

The thickness of the CO layer, which had not been pre- 
viously measured as a function of galactic radius, is roughly 
constant in the region of Wie molecular ring, but increases 
by a factor of about two from R = 7 to 8 kpc. However, in 
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view of -the large error bars associated with the measurement 
for R > 7 kpc, the increase is not yet firmly established. 

An average z, , determined by weighting the thickness values 
for each ring by the area of the ring^ yields <z, > = 59 pc. 

Data on the average thicknesses of other Population I 
objects, which have already been presented in the Introduc- 
tion and are sixitimarized in- Table V.l, show that, with the 
exception of HI and pulsars, the thicknesses of all the Po- 
pulation I objects agree rather closely. Pulsars, v^hich have 
an anomalously broad distribution, presumably as a result of 
their large velocities at formation, will not be considered 
further. 

Only for HI, SNR's, and now CO is detailed information 
on the radial variation of the thickness available over a 
large fraction of the Galaxy. A comparison of these results 
(Figure V-2 ) again shows a general agreement between SNR's 

and CO, while the HI curve has a similar shape but about 
double the value. Thus there appears to be a natural divi- 
sion of the Population I component into two classes; 1) ob- 
jects more or less directly associated with star formation 
or young stars and 2) HI gas and clouds. 
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TABLE V.l 


Thickness of Population I Disk 


Tracer 

Thickness 

Reference 


(z, in pc) 


HI gas 

140 

Baker and Burton (1975) 

HI clouds 

95 

Baker and Burton (1975) 

HII gas 

40-80 

Lockman (1976); Gordon et al (1972) 

HII regions 

40 

Burton et al (1975) 

SNR 

67 

Ilovaisky and Lequeux (1972) 

Pulsars 

330 

Seiradakis (1977) 

CO 

59 

present work 

y-rays 

?? 




Figure V-2. Half -thickness of the galactic disk for CO, 
HI, and SNR's. The HI thickness, which includes both 


clouds and diffuse gas, has been measured at the tangent 
point by Jackson and Kellman (1974). The z, for SNR's 
has been derived from the value of <z > given by 
Ilovaisky and Lequeux C1972) . 
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As was discussed in the Introduction, a similar corre- 
lation of the radial distributions of all Population I ob- 
jects — again except for HI — has already been noted by 
other authors. Based on the radial distribution and the 
limited thickness data available to him, Stecker (1976) has 
proposed that atomic hydrogen be regarded as a new, distinct. 
Population 0. The detailed thickness data presented here 
lends additional strong support to this classification scheme. 

Combining the central emissivity and thickness data, we 

can now calculate the CO surface emissivity of the Galaxy as 

z, e (ir/ln 2} , For R>8 kpc, where there is no thickness 

^ 0 

measurement, we have adopted = 100 pc. The result {Fig. 

V-3) implies that, because of the increasing thickness 

at large R, the molecular ring is not as sharply peaked as 
previous surveys had implied, but extends from 5 to 8 kpc. 

Offset of the CO disk from the con ventional plane 

The position of the CO disk, 2 ^ (R) , is about 40 pc below 
the conventional galactic plane in the region of the molecu- 
lar ring. It was previously noted that the contour plots 
(Fig. IIj-- 2 the gross average of the data over all longitudes 
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and velocities (Fig. IV-1 } , and several munerical tests 
all suggest that such a displacement exists over a v/ide range 
of longitudes, demonstrating that this effect is a general gal- 
actic phenomenon rather than the result of a single large 
complex south of the galactic equator. 

A displacement of the plane of this magnitude is consis- 
tent with the position of the HI plane in the region surveyed. 
Examination of a typical 21 cm latitude-velocity contour plot 
(Figure V-4 ) reveals that the position 

of the emission centroid with respect to b = O'’ fluctuates 
on the scale of the CO plane depression (roughly 
Quigora (1974) and Henderson (1967) have actually measured 
the position of the HI layer and have found that it varies 
about ±50 pc inside the solar radius. 

There has also been a rough determination of the position 
of the SNR layer made by Iloviasky and Lequeux (1972) . When 
the SNR and CO results are plotted together (Fig. V~5 ) , a 

rough correlation is apparent. Lockman (1977) has determined 
that the HII regions also generally lie at negative z. So 
once again it appears that CO acts as a good tracer of the 
early Population I objects. 




Figure V-'4. Distribution of HI emission at il = 25" as a function 
of V __ and galactic latitude, from Kerr (1969) . Note that the 

LjoH 

latitude of the galactic "plane" fluctuates by about 0.3®. 
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Figure V-5, Displacement of the galactic disk from 
the conventional galactic plane as determined using 
CO (present work) and supernova remnants (Ilovaisky 
and Lequeux, 1972) . 
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Comparison wxth Y~^^y data 

Finally, we consider the y-ray data and its relationship 
to CO. As we have seen, the similarity of the y~rcLy and CO 
radial distributions implies that most of the interstellar 
gas is molecular, but, because of the uncertainty of the 
cosmic-ray distribution, the argument is not conclusive. 
However, if, as seems likely, the cosmic-ray disk is substan- 
tially thicker than the gaseous disk, then the thickness of 
the y-ray disk would be an excellent measure of the thickness 
of the interstellar mass . Because the atomic disk is twice 
as thick as the molecular, the y-ray thickness observations, 
when they become available, will give a good indication of 
the relative molecular and atomic masses. 

Summary 

To sum up, a new out-of-plane survey of CO line emission 
in the first quadrant of the Galaxy has confirmed the findings 
of previous in-plane surveys that CO is concentrated in a 
molecular ring about 6 kpc from the galactic center, and has 
furnished determinations of two other important properties 
of the molecular disk as functions of galactocentric dis- 
tance, namely the thickness of the disk, and its displace- 
ment with respect to the currently adopted galactic equator. 
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Comparison with other constituents of the Galaxy has 
confirmed that C0 can- he used as a tracer of young Populra— 
tions I objects and is, in fact, because of its ubiquity, 
the best available tool for observations of this component. 
On the other hand, HI has a different distribution, both in 
2 and in radius, but its thickness does increase at roughly 
the same galactic radius as that of CO. 

This thesis has reported on the first phase of an on- 
going systematic survey of the galactic disk. As additional 
data becomes available, it should be possible to relax some 
of the constraints of the model, particularly azimuthal sym- 
metry, and produce a better picture of the distribution of 
molecular matter. Also, more detailed studies of individual 
molecular complexes, combined with higher resolution studies 
of individual distant clouds done v/ith larger antennas, will 
elucidate the nature of the molecular regions. • 
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APPENDIX A 
Telescope Pointing 


Radio telescopes/ unlike most optical telescopes, normal- 
ly operate v^ithout constant monitoring of their position with 
respect to known astronomical objects. An accurate and depen- 
dable method for determining the position of the telescope is 
therefore essential. In the case of the 4-foot telescope, 
the pointing is further complicated by the absence of astro- 
nomical point sources with sufficient flux at the operating 
frequency to be detected with a 4-foot telescope, and by a 
desire for a high degree of automation of the telescope drive. 

These problems were solved by using accurate computer 
readable position encoders on the telescope axes, calibrating 
the pointing using an optical telescope attached to the ra- 
dio dish, and assuring the colinearity of the radio and opti- 
cal axes. In addition, frequent direct radio frequency checks 
were made against the position of the sun. 

Shaft encoders : 

The position of the telescop'e is determined by Baldwin 

15 bit optical shaft encoders mounted on each axis. These 

15 

encoders give a resolution of 1/2 revolution or 39.6 arc 
seconds (referred to below as one bit) . The encoder accuracy 
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was assured by two testing procedures, one sensitive to 
high angular frequency errors, and the other to low frequency 
errors. - - 

The low frequency errors, which would produce systematic 
pointing errors as a function of sky position, were measured 
by comparing the encoders with the precisely machined worm 
wheels that were part of the telescope drive. The encoders 
were mounted on the telescope (with the antenna dish and 
front end removed) , The axis to be tested was then rotated 
under computer control in approximately 1/1000 revolution 
units using the stepping motor that is part of the normal 
drive system. At each position the shaft encoder was read 
by the computer and compared with the reading that would be 
expected if the worm wheel and encoder were both perfect. 
Figures A-l and A-2 show the difference between the actual 
and expected readings as a function of rotation angle for a 
complete revolution. Note that because 1/1000 revolution 
corresponds to about one half revolution of the stepping 
motors (which have negligible cumulative errors over many 
, revolutions) inaccuracies in the position of the motor steps 
would appear as a very high frequency component in the 
results. 

Figure A-1 shows the results of the test of the 
altitude encoder. The start of both graphs corresponds to 
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ELEVATION ENCODER 



ANGLE 


Figure A-1. Top: Deviation of elevation shaft encoder 
reading from expected value vs position of mount. 
Bottom: Same, except encoder rotated 180° with respect 


to the mount 




AZIMUTH ENCODER 



Figure A-2 . Top: Deviation of azimuth shaft encoder 
reading from expected value vs position of mount. 
Bottom: Same, except encoder rotated 180° with 
respect to the mount. 
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approximately the same position of the telescope mount, but 
different positions of the encoder, which was rotated 180“ 
with respect to the mount between these tests. Most of the 
structure in the plots is the same before and after the 
encoder rotation and is therefore associated v;ith the gear 
on the mount rather than with the encoder. 

Figure A-2 shows the results of tests on the azi- 
muth axis, with the encoder rotated 180“ between tests- The 
source of the shape on these plots could not be established. 
On both axes the peak to peak error is 5 bits of which one 
bit IS due to roundoff effects. The system error is thus 
less than ±2 bits. This error is a combined error due to 
encoder and mount and thus is an upper bound on the encoder 
accuracy. 

To test for high frequency errors, the encoder was se- 
curely mounted in a test jig and the shaft was rotated at 
approximately 1 RPM by a DC motor powered by a regulated sup- 
ply. The encoder was then read at equal time intervals of 
a few milliseconds by the computer and the number of times 
each encoder value occurred was tallied. In an ideal system 
each value would occur the same number of times; deviations 
from this ideal are the result of some bit patterns occurring 
over larger angular widths than others. Figure A-3 shows the 
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ALTITUDE 




2"’ 2'^i 


AZIMUTH 




2^ 2 ‘’-I 

ENCODER READING 


Figure A-3. Test of shaft encoders. Number of occurances of each 
possible encoder reading during test. Ideally, each reading would have 
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occurred the same number of times. 







77 


number of times each value occurred. Histograms of the 
number of values that occurred with a given frequency are 
shown in figure a- 4. From these figures it is clear that 
the deviation in the angular sizes corresponding to indivi- 
dual encoder values is about 1/3 bit. 

Pointing model and initial determination of 
the properties of the mount ; 

The azimuth axis was first set to the vertical to 
within 0.5 arc minutes with a precision bubble level. Then, 
using the optical telescope mounted on the radio dish, an 
initial determination of the properties of the mount was 
made. Following Chauvenet (1960), the pointing of the op- 
tical telescope was modeled by the following equations: 


A = A’ + AA + b cot(z’ + Az) + c csc(z* + Az) 


+ 180° 0(-z' - Az) 

(A.l) 

z = 1 z ' + Az 1 - R 

(A. 2) 

= (j) ’ + A^ 

(A. 3) 

1 = X' + AX 

{A. 4) 


where A = true azimuth 

A’ = azimuth encoder reading 

b = correction due to non-perpendicularity of altitude 
and azimuth axes 

c = correction due to non-perpendicularity of altitude 




0 

NUMBER OF OCCURANCES 


Figure A-4, Data of Figure A-3 histogrammed by number 
of occurances. Ideally, each encoder reading would 
have occurred the Scune number of times and the histo- 
gram would be zero except at a single poxnt. 


79 


axis and optical axis 
z = true zenith angle = 90° - altitude 
z' = zenith angle encoder reading 

({) = apparent latitude (to be used in right ascension 

declination to hour angle, altitude transformation) 

= true latitude of observatory = 40.8101° 

X = apparent longitude 

X'-= true longitude of observatory = 73.9617° 

R = 1 arcmin tan z = correction for refraction 
0(x) “ 1 for X > 0, 0 for x < 0 

and AA, Az, A(f), and AX are offsets. 

^ and X depend on the true position of the vertical (azimuth) 
axis, and AA and Az depend on the mounting of the encoders. 
The absolute value in the expression for z and the 0-function 
are required because two different telescope positions 
correspond to the same sky position. 

Approximate values of AA and Az (the only large correc- 
tion terms) are obtained by sighting on a point (the top of 
the Palisades toward Alpine Tower) v^hose coordinates rela- 
tive to the telescope had been determined from a map (A = 
10,76°', 90-2 = 0.17°) . 

The value of c can be determined by sighting any posi- 
tion on the horizon (z = 90°) in both positions of the tele- 
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scope. In the first case equation A.l reduces to 
A = A’ + AA + c 
and in other positions we get 
A = A" + 180° + AA - c 

where A* and A" are the apparent positions in the two 
positions. Solving for c we get 

c = (A" - A' + 180°)/2. 

Similarly, by observing the pole star in both positions 
of the mount, all the other parameters except X can be de- 
termined, Denoting the azimuth encoder reading, zenith angle 
and time of measurement in the first position by A', z' and t’ 
in the second position by A", z" , and t" equation A.l give to 
first order in b, c, and R. 

Aj^ = A' + AA t b cot z + c csc'z 

A^ = A" + AA - b cot z - c CSC 2 + 180° 

where A^^ and A^ the true azimuths of the pole star at times 
t' and t" as computed from the American Ephemeris and Nau- 
tical Almanac (AENA) pole star table- Approximating z=90°-$> 
and solving for b and AA we get 

b = ^ - c CSC (J) + ^ (A" - A' + 180°)cot(j) 

AA- =-%4A^„ +- A^ - A' A" - 180°)- 
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Using tho constants ^ Sq" (fot times t' and t"), a^, and 
given in the AENA pole star table, the true zenith angle 
at each position can be expressed as 

z' + Az - R = (1) “ a ’ - a, ~ a„ 

^ o 12 

\z" + Az[ “ R= “Z" - Az - R 
“ ■ ®o“ ■ " ®2 

Solving for z and (j) we get 

Az= ^[(“Z") - z' + a^*' - Sq ' ] 

= ^[(- 2 ") + z' + a"+a '] +a-i +a„-R 

o o X 2 

As a result of the procedures described above, the 
values of (ji , AA, and Az were determined to sufficient accu- 
racy that it was possible to easily find stars in the optical 
telescope. Since the value of b, the non-perpendicularity 
of the mount axes, was less than 1' of arc, it was assumed 
that it could be ignored. This hypothesis was later veri- 
fied more accurately by the star pointing measurements. 

Alignment of optical and radio axes : 

Using radio measurements, the only position that can 
be easily loca.ted with the 4-foot telescope to high accuracy 
is the center of the sun. Because the limb of the sun is 
extremely sharp at millimeter wavelengths, it was easy to 
place the radio axis of the telescope on i^he center of the 
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sun by measuring the position of several points on the 
lunb. With the telescope tracking the center of the sun 
on its radio axis, the axis of the optical telescope was 
adjusted until its crosshairs also lay at the center of the 
sun. This was done to a fraction of an arc minute by fo- 
cusing the solar image on a circular "target" (a piece of 
polar coordinate graph paper will do) the same size as the 
solar image and adjusting the optical telescope until the 
sun was centered on the target and the image of the cross- 
hairs was at the center of the target. 

Star pointing ; 

With the radio and optical axes aligned and approxi- 
mate parameters determined, the model of equations A . 1-A , 4 with 
b = c = 0 was used to point the telescope, under computer 
control, at a number of stars scattered throughout the sky. 
After a star was located in the optical telescope, the 
telescope was moved until the star was exactly centered 
on the optical crosshairs and the apparent star position 
was recorded. By varying the parameters of the model to 
minimize the least-squares deviation between the coordinates 
given by the model and the true positions of the stars, im- 
proved model parameters were found. It v/as determined that 
b and c could be taken to be zero without introducing any 
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significant error. 

With the improved values of AA, Az, AX and A4» inserted 
into the model, the deviation of each star was again computed 
and the RMS deviation was calculated. The final test included 
31 stars. The maximuin error at any star was 2.2' and the RMS 
error was 1.1'. Figures II-2 and II-3 show the resi- 

dual errors as functions of azimuth and elevation. 

Solar pointing ; 

Additional frequent checks of the pointing •were made 
against the position of the center of the sun using radio 
measurements. An automauic computer program scanned the 
telescope across the sun rn azimuth and elevation, computed 
the posxtion of the center of the sun by finding the 
points of maximum symmetry of the intensities measured 
along each axis, and computed the error from the expected 
(true) solar position. These checks gave peak errors of 
2' and typical errors of about 1.5'. 
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APPENDIX B 
Antenna Pattern 

The antenna was designed to be diffraction limited . Be- 
fore delivery the manufacturer measured the surface of the 
main reflector and found a peak error of 0.0014 inch. This 
corresponds to X/70 at 2.6 mm, much better than the usual 
engineering standard of X/40 for diffraction limited systems. 

In order to verify that the dish actually met this spe- 
cification and also to determine the tolerances permitted 
in positioning the feed horn, a series of calculations were 
made and were compared with actual measurements of the an- 
tenna pattern. The telescope was tested at 115 GH 2 : against a 
transmitter across the Columbia campus at a distance of 420 
meters corresponding to .72 D^/X. In antenna testing, the 
transmitter is conventionally placed at a distance at least 
2 d /A from the antenna; however, Johnson et al. have shown 
that patterns very similar to the far field patterns can be 
obtained at as little as D^/8A if the telescope is refocused 
on the near field transmitter (Johnson et al. 1973) . 

The theoretical antenna patterns were calculated using 
scalar Kirchhoff diffraction theory. All the calculations 
described here are for the "equivalent parabola" (Hannan 
1961-)-,— the -locus of the i-ntersee-tions of the incomi-ng— rays 
with the extensions of the corresponding rays reflected from 



85 


the secondary (see figure B~l) . This equivalent parabola is 
the same diameter as the actual primary reflector and has 
the same focal length as the Cassegrain system. Near the 
focus, the pattern and focusing properties of this single 
parabola are very close to those of the original Cassegrain 
parabola and secondary. 


The field of an arbitrary aperture in the Fraunhoffer 
region is given by (Silver 1964) : 


U{0,(j)) = (i/XR)e g(e,(})) 
v/here g ( 0 , (J)) = "^aperture 
k = k sin 6 cos (|) 
ky = k sin 0 sin (p 
k = 2 tt/X 


(k X + k y) 
X y-^ 


dxdy 


and F(x,y) is a complex function which gives the field over 
the aperture . 

F(x,y) is approximated in this calculation by zero 
outside the dish and on the blocked portion of the aperture 
and elsewhere by’'' 

F(x,y) = F(p) = (1 - (p/a)^)x (1 - E) + E 
where E is the ratio of the fD.eld at the edge of the dish 


*soe section on "Horn Pattern 
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Figure B-1. Definition of equivalent parabola of focal 
length £ and radius a, p is the radial coordinate in 
the plane of the equivalent parabola* 





87 


to the field at the center. Then 
g( 0 ,tf>) = - Sg ( 0 ,<j)) 


where 


P F(p)e 


ikpsinOcos ( 4> - <() ' ) 


is the field of the unblocked aperture and g_ has the same 
integrand but is integrated over the blockage. The blockage 
is approximated by a central circle representing the secon- 
dary and three rectangles representing the arms (see figure B-2). 
Thus is given by g^ = g^ + g^ v;here g^^ is an integral over 
the central blockage and g^ is an integral over the arms. 


Using these approximations the integrals can all be done 
analytically. Letting r = p/a, Ujj = irx sin6/A and 
f (r) = F (ra) we get 


g. = 27 Ta® 

tJi i-J ZB, 
U 2 .a 

^ p + J?„(U 5 .a 1 
U|a ' 

9 c ( e,<p]=2TTh^ 

'jJUabl 


Ua.b 


i-E 




Jkf Uab 1 





wi-i 


J ^ 


ELW sinCUw 

Ut,cosC^-d:) Uw sin -cc) 


X exp (i 2 Ut,+u/ 2 . cos( 0 -cc)) 

+(l- E) 1 _r*) expri2'-Ja'^°sC0-oc)rldr 

' ' Uw Sin(^-cC) / It ^ 





b = 7.6 cm 
a = 6 1 .0 cm 
w= 2.5 cm 
L=53.3 cm 

Figure B-2. Aperture and bloclcage used in calculation 
of antenna pattern. The circle of radius a is the 
aperture and the sh ided areas are blockages. x, y, 
and p define the coordinate system used in the calcu- 
lation . 
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The pattern intensities are proportional to |gj . 
Figure B-3 shows the intensities in the altitude and 
azimuth planes normalized to O dB at the central maximum 
for 10 dB edge illumination. These graphs also 

show the pattern that would be obtained with an unblocked 
circular aperture and with only the central blockage. Note 
that without the blockage the pattern is azimuthally symme- 
tric in this approximation. Figures B-4, -5, - 6 , and -7 
show several sets of measured points along v/ith the calcu- 
lated pattern. 

Horn Pattern 


Figure B-8 shows the dimensions of the two available 
horns. Since the phase errors across the mouths of these 
horns are small (<A/6 ) , the horns can be modeled to good 
precision by a rectangular waveguide antenna. The field in 
the central portion of the pattern is (Silver, 1964) : 


E, 


Sin u. 


u,. 


Se<%> 


E, 


TT 


cos u 


2 2 
u - ir /4 
a 


= 9h 


where u = sin 0 
X A 



Cl) 



Relative Pov,-er (db) 



Pattern Angte (Minutes of Arc) 


Figure B-3. Calculated antenna pattern with 10 dB edge' illumination. 


VQ 

O 




Pa)iem Angle (Minutes of Arc) 

Figure B-4 . Antenna pattern in horizontal plane with 10 dB edge 
illumination. Solid lines are the calculated patterns. Dotted 
lines are the measurements. 
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1 


/ ’'A 


Pattern Angle (Mhetes of Arc) 

Figure B-5, Antenna pattern in vertical plane with 10 dB edge ilium 
ination. Solid lines are the calculated patterns. Dotted lines are 
the measurements. 
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(0 db HORN 

13 db HORN 

A 

0.618 ±0.003 

0.7!5± 0.003 

B 

0.431 ± 0.003 

0.486 ± 0.003 

L 

3.50 ±0.03 

3.50 ± 0.03 


Figure B-6./ Dimensions of the tv/o available feed horns. 
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and En E, are the E-fields in the E and H planes 

0 (p 

of the waveguide. 

For computational purposes the field was approximated by 


g(p) = E + (1-E) X [1 - {p/a) ] 


where E is the field at the edge of the main dish. Figirre B-7 
shows the aperture illumination functions g^, and g for 
the 10 dB horn. All these cruves lie within 1.5 dB of each 
other. 


Feed Horn Positioning 
A.) On~axis defccasing. 

The predominant effect is a phase error across the 
aperture given by 

'P = 2tt i (1 - cos6) 

where d is the amount of defocusing and 0 is the aperture angle in the equiv- 
alent parabola (see Fia. B-1) . For a long focal length system 0«1, thus 

^ 2ird 0^ 

^ " A 2 



If we limit to Ymax < tt /4 , the effect of the phase error 
wi-ll—be— smal-l-i Thus -we must have 

d X = 0.78 inches 




01 2 345 6789* 10 


Figure B-7. Aperture illumir.ation functions for the 
10 dB horn in the E- and II~planes and the approximate 
function used in the calculations . 
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Experir.icnt confirms that it is adequate to focus to ± 3/4 
inch (see Figure B“10) , 

B.) Off^-axis defocusing 

The result of moving the horn off axis is well covered 
in the literature (Ruze, ] 965; Sandler, 1960; Kelleher, 1952; 
Lo, 1960) . The major effect is a shift in the position of 
the primary lobe of the antenna system. If, in the equiva- 
lent single reflector antenna, the horn is shifted an angle 
a with respect to the center of the parabola, the beam will 
shi£t by 3 = Bet where B is known as the "beam deviation 
factor" and is a function of f/D. For an f/2.8 system, we 
can take B = 1 to high accuracy. 

In addition, when the offset becomes large, there is 

an asymmetric coma distortion of the -antenna pattern. Ruze 

(1965) shows that there is negligible (<1 dB) change in 

2 

If the horn is shifted by a 
distance d (d<<f) then ct~ d/f. The criterion is then 
d<2A (f/D) ~ 4,2 inches. Since the feed is easily positioned 

to this accuracy, no distortions of the pattern should be 


the pattern shape if a<~2[~ 


seen . 




Figure B-8. Focus characteristics with the 10 dB horn. The solid line 
shows the antenna pattern with the horn at the focus. The two dotted 
lines show the pattern with the horn about one inch either side of the focus. 


VO 
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C. ) Misalignment of the Feed Horn 

I-f the feed horn i-s m-r-sa-l-i-gned-/ such that the max-i-mum- 
of the horn illumination is no longer in the center of the 
aperture, the phase will still be constant across the aper- 
ture. If the maximum of the illumination falls at the point 
(0, pa) in the aperture plane (0_<p£l) , the approximate il- 
lumination function becomes 

F(x,y) = E+(l-E)x[l-(x^+(y-pa) 2)/a2] 

Antenna patterns have been calculated using this illumina- 
tion function for the dish with central blockage (no arms) . 
Figure B-9 shows the results for various values of the horn 
misalignment angle a. The graph shows that for a = 1°.0 
the effect is small, but it grows rapidly for a > 1°.0 . 

This is consistent with our experimental findings. 

Accuracy of Dish 

Ruse (1966) reviewed the effects on the pattern of 
random errors in the figure of the reflector. He shov7ed 
that for small errors the central part of the pattern is 
given approximately by the power sum of the no-error pat- 
tern and a broad pattern of scattered radiation. The ratio 
of' the scattered radTafibn to the no error orf axis maximum 

1 2 c 2 

radiation is R = -- (-g-) y- v;here n is the aperture effi- 
ciency, D IS the diameter of the 


aperture, <5 is the error 



Relative Power (db) 



Pattern Angle (Minutes of Arc) 


Figure B~9. Calculation of the effect of horn misalignment 



xoo 


in the figure of the dish expressed in radians at the opera- 
ting wavelength, and c is the typical correlation length of 
6. Because the error pattern is very broad compared to the 
mam beam, the first null in the pattern relative to the on 
axis maximum must be greater than R. If the measured inten- 
sity at the first null is e relative to the maximum, then 

e>(l/n) (2 c/D)2 $2 ^ (2c/D)2 6^ 

Examination of typical patterns shows that e < 10“^*^. Thus 

Sms E V 6'2 < 2 . 10”^ 5 . ~ i (inches) . 

~ c c 

There may be errors of a significant fraction of a V7ave- 
length for c - 1 inch, but for large c, errors must be quite 
small. For example, for c = 10 inches, 6rms < 1/10 radian 
= 1/63 wavelength = 0.0016 inch. "For a correlation length 
equal to the diameter, we get 6rms < 0.0003 inch. This 'is 
consistent with the measured peak error of 0.0014 inch. 

Conclusion 

The antenna pattern has been calculated and compared 
with actual measurements. It was determined that both the 
antenna surface and the feed horn alignment were adequate 
to produce a diffraction limited antenna pattern that agreed 
well with calculation. The bearawidth was S arcmin full 
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v/idth at half iriaximum and the sa delobes are more than 15 dB 
below the central maximum. 



APPENDIX C 


FITTING THE MODEL TO THE DATA 


The clumpy nature of the CO spectra strongly suggests 
that CO IS concentrated in well defined molecular clouds. 

We have not attempted to predict the positions of these 
clouds in detail, but rather, have assumed, as was ex- 
plained in Chapter IV, that the clouds are randomly dis- 
tributed in a vfay that can be described by an effective 
emissivity, s, which is a function of galactic radius and 
distance above the galactic plane. As was shown, once e 
is given, it is possible to calculate the expectation value 
of the antenna temperature, T*, as a function of galactic 
coordinates and velocity. (This is worked out in more 
detail in Appendix D, including the effects of beam 
smoothing and the finite resolution of the spectrometer.) 

This appendix shoves how we can select a model function 
for e which, in a sense to be described, predicts antenna 
temperatures which best fit the actual observed temperatures. 
In addition, we will discuss the accuracy of the fitted 
parameters . 
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DGrivatiori of least-squares fitting procedure 

We begin by follovjing the usual procedure for fitting 
a model function to experimental data. First, using our 
knowledge of the physical situatxon, we choose a general 
form for e(R,z) which depends on several parameters — for 
example, and in the model of Equation IV -8 — 

whose values are to be determined by the fitting. 

From the model effective emissivity, we can bhen cal- 
culate, as a function of velocity and galactic coordinates, 
the expectation value of the antenna temperature and its 
dependence on the parameters. Explicitly, we write 
T* = T* (X,A) , where abbreviate the coordinates and 
velocity by X = (ii,b,v) and the entire set of parameters 
by A = ' ^2'’’*' ^ e.g. in equation IV -8 we have 

A = (e_, z, , z ) and N = 3 . 

0^0 

The usual practice would now be to find the "best" 

value of the parameters by the least-squares procedure; by 

varying the parameters. A, so as to minimize the expression 
s _ 2 

ill (^i-f A) - T*^) , where the experimental data 

are 'a series of s measurements of the antenna temperature, 
T*^ (1=1,..., s), made at the positions and velocities 
^i ~ ^^i' ^1' However, this method really applies 

only if the statistical fluctuations in T*^ have a 
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gaussian distribution. In the case at hand, the distri“ 
button of the fluctuations is unJcnown and may be very 
different from a gaussian. For example, to take a very 
simple case, assume that the Galaxy is filled with identical 
clouds, each with the same radiation temperature, . The 
number of clouds, n^, observed in a given spectrometer 
channel is then Poisson distributed. If we further assume 
that each observed cloud has the same beam filling factor, 
f, then the antenna temperature in the channel is given 
by T*=n T f. Ifn <1 then T* is unlikely to have 
any value other than 0 , f T^ , or 2 f T^ . If we con- 
sider the variation in antenna temperature with, for 
example, galactic longitude, then, in contrast to T*{X,A), 
which is smoothly varying with longitude, the actual 
observed is a quantized function which jumps between 

0 and f . (See Figure C-1.) It is not apparent how 
we should fit the smooth T'^ to the very different looking 

iji* 

i 

In fact, V76 do not know whether the fluctuations in 
T* are gaussian, poisson as in the example, or perhaps some 
intermediate distribution. However, we will shovj that the 
fitting procedure is similar for both the gaussian and 
poisson cases and that the assumption of either form 
leads to similar results for the fitted parameters. 



Antenna Temperature 
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Figure C-1, Comparison of model and observed antenna 
temperatures for the simplified model described in 


the text 







3 OS 

Follovjing Janossy (1965) , the derivation of the 
procedure for least-squares fitting will be summarized 
briefly for both the poisson and gaussian cases. The 
theory of least^squares- fitting is usually derilved/'. 
from the maximum likelihood principle. The probability of 
obtaining the measurement results Y = / Yg) of a se- 

ries of s measurements made at the points (independent vari- 
ables) X = (X^, .../ is given by P(Y?X,A), where A stands 

for the entire set of N parameters that are to be determined, 
i.e., A = (a^, ...f - The best choice for A can be taken 

to be the value that maximizes P or, equivalently, the like- 
lihood function, P. If the y^ are independent random va- 
riables, we can write 

^ A) = S Cfip(5^.; X., A) (C-3) 

i 


where p gives the probability density for a single measure- 
ment. The best choice for A is then given by the solution of 
the equations 


AfeiJP 
S a 


: = I 


9^ P(yj5 , A) 


d a 


= 0 


T = 1, . 


J1 


(C-4) 


The variance of a can be shown to be 
i 


-X 

var(a ) = G ^ where G 

tt Xt 


3 0n P 

a a, S a 
X T 


(C-G) 
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If the measurements have a Gaussian distribution^ 


P(y ; X , A) = - 


3 r / ^ -j-j 2 \ “1 

~ exp [- (y. -y^) / (2 o - )] 


(C-G) 


where y. = ^ is the expectation value of jn . Then, 


0n P 


- I I 


fri - y^r 


+ constants 


(C-7) 


paid the juaxt^iiizaiion of j®?i P is immechalcly seen to lie the usual least- 
sqimrcs condition. Writing the cejuations explicitly, 


P 
d a 

T 




0 


(C-8) 


G 


tX 


em P 

B a,, B a 

T X 


V 


2 

^i 


B a„ B a 


(C-9) 


Por the Poisson disLiibution, 


~ Y y- 

' I ~ -^1 


p(V.; X, A) - e y. / y. 


(C-10) 
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TJicn 


SmV - Z (y. y- " y- ) '*■ containing A) , (C-11) 

i ^ 


The miixirmim likclilioocl conditions are 


I 





(C-12) 


G 


tX 





(C-13) 


These egtiations are similar to Eq, (C-8) and Eq. (C-9) ; 
the only difference is that is replaced with y^. For 
poisson statistics, y^ is in fact equal to the expected 

variance of y . . Hov/ever, y is a function of A to be 
determined from the fit, whereas in the gaussian case 
is a constant which must be known before the fit is made. 

The two sets of equations are thus actually quite different 
in spite of their apparent similarity. 
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As outlined, by Janosky (p. 258) Eq . (C-12) can be 

solved to adequate precision by a minor modification of 
the usual least squares formalism. (2\ctually^ Janossy 
presents a more general case.) Initial values, , are 

estimated for the parameters, A. Then Eq. (C-8) , the 
normal least squares equation, is solved with held 

fixed at the value y(X^, From the solution new 

values, , are obtained and the process is repeated 

with now set equal to y(X^, '). This process is 

repeated until it converges. 

In the usual terminology of least-squares fitting, 

t h 

l/a^^ is called the weight of the i^ poinu; the points 
with the smallest uncertainties are given the largest 
weights. Our iterative procedure thus corresponds to a 
series of least-squares fits with the V7eights adjusted 
at each sbep to correspond to the best estimate of the 
uncertainties as determined from the previous step. That 
this procedure works is intuitively reasonable since the 
results of a least-squares fitting generally do not 
depend strongly on the vzeighting used. 

Our fitting algorithm was checked using a stocliastic 
simulation of the Galaxy simlini to the one described in 
Chapter IV. The only differences were that the width of 
the “molecular rjng" was doubled by setting = 2 kpc 
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(seo Equation IV'” 9 ; and the number of clouds was 

decreased so that at the peak of the CO distribution there 
was only about one cloud per resolution element of the 
simulated spectra. The resulting spectra were fed into 
two versions of the least-squares fitting procedure; 1) the 
standard code with constant weighting — correct for gauss- 
ian statistics, and 2) the modified code — correct for 
poisson statistics. As can be clearly seen in Figure C- 
the poisson case produced considerably improved results. 

The error bars are a much better representation of the 
actual errors, and we \?ere able to determine the para- 
meters for two radii at which it was not possible using 
gaussian statistics (R > 7.5 kpc) . 

In practice, the procedure we have outlined must be 
modified slightly to take account of additional fluctu- 
ations due to receiver noise. To illustrate the problem 
we return to the example of many clouds, each with the 
same beam filling factor, f, and temperature, . When 
we observe a point far from the galactic plane, the number 

of clouds in the beam, n , will be small. Thus the 

c 

variance, , of n , which is equal to n for poisson 

c c 

statistics, will also be small, and the weight of the point 
will be larg^_ If . th,e__rejo_e.iv.er._ v;.er_e„no.is.e rfr.ee,, „ this, would 
be correct, for the detection of emission far from the 
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R(kpc) 

Figure C-2. Results of fitting to a simulated galaxy 
with a small number of clouds in each resolution element. 
The results on the left were obtained using the standard 
least-square code with constant weighting? those on the 
right using the modified code for poisson statistics. 

The dashed lines { ) indicate the underlying distribu- 

tion used in sunulating the galaxy. The dash-dot lines 

(-• •-) show the original parameter estimates used in 

the fatting. 









plane would be an extremely unlikely event. In reality, 
once the fluctuations in antenna temperature due to the 
cloud statististics fall below those duo to receiver noise 
increasingly larger weights should not be assigned. 


To take this into account, at each step of the itera- 
tion the weights were chosen to reflect the best estimate 
of the fluctuations, including both receiver noise and 
cloud statistics; that is, the variance in the antenna 
temperature was taken to be 


o 


i 


2 



+ CT 


C/1 


(C-14) 


where is the reciever noise, and a . ^ is the variance 
N C , 1 

due to cloud statistics. Using our simple model again, we 
estimate a ^ = (f T ) ^ n , , where f T is constant and 
for poisson statistics, n is the variance in n . . If 

we estimate the contribution of a single cloud to the 
antenna temperature to be f = 1 K, then 


a 


c 




n 

c,i 


= 1 K (f T n . ) 

c c,i 

= (1 X) T* {X^, ) 


(C-15) 


where we have used the best estimate of from the previ- 
ous iteration. So, putting the results together for the 
poisson case, we used the v/eightinvj function 
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1/0^2 = [ ^ (2. K) T*(x^r . (C-16) 

In the limit of = 0 , this is just the exact result 
obtained in Equation C-^12 . 

We have actually fit the data both using the weight- 
ing as just described and also using a constant weighting 
function, which would be correct if the antenna tempera- 
ture fluctuations were caused primarily by receiver noise. 
The results for these two choices, when applied to the 
model described in Equation IV-8 (emrssivity is a 
gaussian in z) , are shown in figure C-3. In the region 
of the molecular ring, both choices of the weighting 
give about the same results. The only significant differ- 
ences appear in z, for R > 6.5 kpc, where the increase 
in thickness is somewhat less for the constant weighting. 

The results given in the main body of this thesis were 
all obtained using constant 'weighting. Although this 
ignores the problem of counting statistics, we feel that, 
in view of the demonstrated insensitivity to the choice 
of weighting fiinction, it is preferable because it does 
not introduce additional parameters, such as f T^, whose 
values are extremely uncertain. 



Figure C-3 . Effect of choice of weighting function, - 
Left; Dispi acement , half thickness and central emissivity 
of CO disk determined with a. ^=T* CX. ) + T^t^* 

Right7' Same parameter^determined vjith~dT^'=: constant. 

See text for details. 
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Acc\;racy of the p ar ameters 

We turn now to a brief descussion of the accuracy to 
which the parciineters can be determined. The variance of 
the fitted parameters could be determined using Equations 
C-9 with the set equal to the values used in the last 

iteration of equation C-8. However, the variances of the 
parameters, unlike the values of the parameters, are sen- 
sitive to the values of the 0 ,^ , which are not accurately 

1 

known. Therefore, we used instead the usual procedure 
for finding "internal errors" when the fluctuations 
in the data are not accurately known. That is, 

2 

the 0 ^ were first scaled by a factor k so as to make the x 

of the data equal to its expectation value: 
si ~ 

2 (T*.- T*(X.,A)) = s-H, These scaled values, 

i=l (ka^)^ ^ "" 

2 

(k 0 ^) , were then used in eq. C-9. 

The preceeding discussion has assumed that the T*^ are 
all uncorrelated. In fact correlations probably exist because 
the linewidth and spacial extent of a cloud cover more than 
one resolution element of the survey. To include these effects 
properly, additional parameters describing the properties of 
individual clouds would have to be introduced. Because the 
quantity of data now available does not seem to justify such 
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complex models, the correlations have been ignored. The ma- 
jor e-ffeet of th-x-s gross approxx-mation will be to underesti- 
mate the errors on the fitted quantities. 

The correlation can be roughly estimated by examining 
the spectra. A typical feature seems to cover about 3 chan- 
nels (7.8 km sec and 2 resolution elements in latitude 
(0.5°) or a total of 3 x 2 = 6 resolution elements. The er- 
rors could, in the worst case, be underestimated by about 
/6 . In practice, however, the errors bars probably reflect 
systematic errors in the model and are relatively insensi- 
tive to the quantity of data. The error bars should thus 
be interpreted as only a rough guide to the size of the ac- 


tual errors. 
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APPENDIX D 

CALCULATION OF MODEL ANTENNA TEMPERATURES 

In order to compare a model emissivity function with 
the survey data, it is first necessary to convert emissivity 
as a function of galactic position into antenna temperature 
as a function of telescope position and frequency. In doing 
this calculation several approximations have been made. These 
will be described jn this appendix. 

If, as we have assumed, the Galaxy is azimuthally sym- 
metric, then the model specific volume emissivity (or emis- 
sion coefficient) of CO can be written e~(R,z) where R is 
galactic radius, z is height above the galactic plane, and 
V IS the frequency as measured by an observer moving with 
the gas at location (R,z). To evaluate the intensity at the 
earth it will be convenient to use the earth-centered spheri- 
cal coordj nates (2.,b,D) where I and b are the usual galactic 
coordinates and D is distance from the earth. R and z can 
be expressed in these coordinates as 

R(S,,b,D) = (d2 + _ 2r^ D cos b cos (D-1) 

2 (b,D) = D sin b 

where = 10 kpc is tiie distance from the earth to the 
o 

galactic center. Because the gas as in general moving with 
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respect to th.e local standard of rest (LSR) ^ it is usedEul to express 
V in terms--©f the ©laserved Doppler shifted-frequency- v: 

V = v(l + V/c) (D-2) 

where V is the line-of-sight velocity of the gas with res- 
pect to the LSR If the Galaxy is in circular rotation 
then (Burton, 1974) 

V = R^[o)(R) - w(R^) sin 1] (D-3) 

where w(R)is the angular velocity of the gas. We have used 
Burton's (1974) formula for w(R): 

R w(R) = 250.0 + 4.05 (10-R) - 1.62 (10-R) ^ (D-4) 

The specific intensity as measured from the earth is 

CO 

I (J!,,b) = / (R,z) dD (D-5) 

V o V 

where v, R, and z are given by equations (D-1) and {D-2). We 
are interested in the intensity averaged over a single chan- 
nel of the spectrometer; 

_ v+Av/2 

= (1/Av) f I , dv' (D-6) 

v-Av/2 ^ 

where Av is the width of the channel centered at v. Substi- 
tuting equation (D-5) into equation (D-6) we get 

_ v+Av/2 

i., = (1/Av) / / c~, (R,z) dv' dD (D-7) 

0 v-Av/2 

where v' is derived from v' by cquarion (D-2). 
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In the absence of any iii formation of the line shape we 
assume that the line is sufficiently narrow compared to Av 
that we can make the approximation 

e^(R,z) = e(R,z) 6(v-v^) (D-8) 

where e(R,z) is the emissivity integrated over the line and 
is the rest frequency of the CO transition. We can now 
evaluate the inner integral of equation (D-7) as 

v+Av/2 e(R,z), v-Av/2<v <v+Av/2 

/ 6;^, (R,z)dv' - { ° (D-9) 

v-Av/2 0 , otherwise 

Substituting into equation (D-7) we get 

= (1/Av) /e [R(£^b,D) ,z (b,D) 1 dD (D-lO) 

where the integral extends over all D satisfying 

1 V ~ v^I < Av/2 (D-11) 

or, substituting for v from equation (D-2) 

i V [1+V(£,R) /c] “ ^qI < Av/2 (D-12) 

where R = R(t,b,D) . 

BV 

As — < 0 everywhere (ah least for the model rotation 
BR 

curVe wc are using) , it is clear that the positions along the 
line-of-sight that satisfy equation (D-12) must lie between 
tx/o radii R and R . (See figure D-1) . If R < R , where 
R^^ = R^ sin £ is the galactic radius of the subcentral point, 
then in the inner part of the Galaxy (i.e., there are 
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two segments along the line-of-sight that satisfy equation 

(D"12) . If then the two segments merge. 

However if we redefine e R in this case, then it reduces 

< sc 

to the previous one. If R < R then there is no solution 

^ sc 

of equation (D-12) along the line-of-sight and = 0- Thus 

we need only consider further the first case, R < R^. 

sc < 


We now evaluate in equation (D-^IO) . The integral 
can be split into two pieces, one covering each of the two 
segments between R^ and R^ . In each segment the function e 
was replaced with its value at the point along the line of 
sight at the radius R = h{B. + R ) . If D, and are the 

9. V y ^ ^ A 

distances to these two points and AD is the identical length 

of each segment then 

2 



(AD/Av) E s[r (£,b,v) 

•I 1 3 


z (b,D^)] . 


(D-13) 


The replacement of R by causes only a small error 

in I^ but results in an enormous simplification of the final 

least-squares fitting. For the spectrometer that was actually 

used, - R^ < 0.4 kpc over the entire area of the survey. 

Thus, in substituting R for R in equation (D-10) , e is eva- 

avg 

luated at a radius less than J| (R^ - R ) =0.2 kpc from where 
it should be. This is less than the resolution of the models 
we have used, and in any case larger errors result from the use 
of a circular rotation model. With this approximation, / 

and ultimately the model antenna temperature, depend on the 
value of e at only a single value of R. Thus, when, in the 
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final models, the Galaxy is split into concentric rings, 
the data points can be grouped according to which ring they 
depend upon. This lack of "mixing" between the rings allov^s 
the parameters in each ring to be least-squares fitted in- 
dependently, thus vastly reducing the amount of computation 
needed. Finally, to get the model antenna temperature v/e 

must convolve I^(Jl,b) with the beam of the antenna and mul— 

2 2 

tiply the result by c / (2kv ) to convert to the usual 
radio astronomy units: 

1 C.2 

T*<)l,b,v) = 5 - G(Z',b'; l’.b)da' (D-14) 

4tt 2 k Av 

where G (£' ,b' ; Jl,b) is the gain of the antenna (relative to an 

isotropic radiator) in the direction (£',b‘) with the axis 

of the telescope pointed in the direction (£,b) . If we de- 
2 2 

fine s* = sc / (2 k v ) , a convenient unit of e for radio 
astronomy, and substitute equation (D-13) into equation (D-14) 
we get 

AD ^ 

T* = E // 6*[r ^^^(a',b',v),z(b',D. {f,b',v))j 

47tAv i=i 4 tt ^ 

' ,b' ;£,b) dfi' (D-15) 

■Because the telescope beam is very small compared to the 

SjC,aLe._jon...whi_ch we will attem pt to measure radia l fl uctuations 

in the emissivity, it is safe to evaluate R and D. only 

avg 1 

~t the center of the beam. 
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The gain function of the telescope can be accurately 
itiodcled by 

G(£' ,b' ;£,b)=^(4/l2)exp[-( (£ -S- ' ) (b-b' ) ^) /A^] CD-16) 

where X is chosen to produce a halfwidth at half maximum of 
8 arc min. Since the beamwidth is much less than one radian 
and b << 1, we can replace the spherical integral by 


// dfi ^- / / 

TT -w -« 


db 


and also make the approximation 

z(b,D) ~ bD. 


(D-17) 


(D-18) 


The integral over 5,' can now be done analytically with 
the result 


•j** = 


I \ 2 2 ' 


exp[“ (b-b' ) ^/X^]dh' . 


CD-19) 


In practice r we have considered t\^o forms for the func- 
tion e* — gaussian and exponential in z. In both cases the 
remaining integral over b’ can be done analytically. In the 
gaussian case we take 


e*(R,z) = e^(R) exp[-(z-z^(R) ) ^/a(R) (D-20) 


where a(R) = Zj^(R)/(ln 2) ^ 
Then 


ADc 


T* = 


2 

S 


a 


Av 


1=^1 (c^ + P^)-^ 


oxp [- 


(z .-z ) ■ 
1 o_ 

a^ + 3 . ^ 
1 


} (D-21) 
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where 6. = AD. 

1 - X 


In the e^ponentxal case we take 


e*(R,z) = exp [- jz“Z^(R)l /a (R) ] (D-22) 


where a (R) = Zj^/(ln 2) 


Then 

AD e 2 ry 

T* = o_ E exp [ 3 . /(4a^)] x 

2AV x-i ^ 


(D-23) 


P . ifi • 

{exp[(z^-z^)/a] erfc ( — ^ + — — ) + exp[- (zj^-z^) /a] 

2a 

3 . z. 

erfc C —i ~) } 

2a 


where 3 . and z . are the same as in eg (D-21) and erfc is the 
complementary error function. 


Equations (D-21) and (D-23) are the actual exprcssxons 


that were used in the computation. 
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